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ABSTRACT 



Fluctuating temperature structures in the atmosphere 
induce phase perturbations in a propagating laser beam. These 
turbulent conditions occur throughout the atmosphere and cause 
the laser beam to spread and alter its centroid. There are 
several methods to measure the parameters of optical 
turbulence in the atmosphere, but few that will determine them 
as a function of altitude at all levels. One method of 
measuring altitude profiles of turbulence is with a 
temperature probe launched via a balloon system. 

This thesis involves the development of a differential 
temperature probe sensor to measure the temperature 
fluctuations at all altitudes in the atmosphere. In addition, 
it investigates the effect of solar heating on the probes in 
the atmosphere and the subsequent effects on the measurements. 
A validation of the probe system was made by a comparison test 
with an acoustic echosounder developed earlier. In addition 
to validating the probe system, the absolute C T 2 analysis of 
the echosounder was confirmed. 
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I . INTRODUCTION 



Atmospheric conditions will cause severe degradation 
along the optical path of a ground to space weapon or 
surveillance laser. There are several causes for this 
degradation, they are 1) absorption and scattering by aerosols 
such as rain and clouds, 2) distortion by thermal blooming 
and 3) distortion by atmospheric turbulence [Ref. 1] . 
Absorption and scattering can be controlled by varying the 
wavelengths of the laser and having multiple sites to insure 
at least one has a cloud free line of sight. Thermal blooming 
is the heating of the medium, through which the laser beam 
propagates. It is due to absorption of the radiation by 
molecules and aerosols and the consequent distortion of the 
beam due to density reductions brought on by the heating. 
Choosing approximate wavelengths of the laser which have low 
atmospheric absorption in the atmosphere reduces this effect. 
Atmospheric turbulence is difficult to deal with since there 
is no way to avoid it. 

The major effect of turbulence on an optical beam is the 
limitation of the mutual coherence lengths. For example, an 
average ground to space coherence length for the atmosphere 
is on the order of 5 cm therefore a ground based laser having 
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a 4 meter diameter mirror will deliver less then 1/1000 of its 
original power onto a target. Adaptive optics provides a means 
for reducing these turbulent effects. It corrects for the 
effects of turbulence by altering the wavefront 
characteristics of a beam [Ref. 2] using deformable mirrors 
or nonlinear optical materials such as Barium Titanate. 

Measurements of the turbulence from the surface to an 
altitude where the turbulence is no longer significant 
(approximately 30 km.) are needed for several reasons. The 
most important is to be able to characterize the turbulent 
profiles of the atmosphere at different locations in order to 
select the best site for a ground based optical or 
surveillance system. Each site will have a different 
turbulence profile since it depends on the local geography as 
well as the upper atmospheric conditions controlled by general 
meteorological patterns. Another reason for the measurements 
is to determine vertical distribution of the turbulence 
parameters, such as the coherence length, which affect the 
design of adaptive optics systems, or the signal processing 
transformations in a surveillance system. 

There are several instruments for measuring turbulence 
parameters, some of which measure the index of refraction 
structure parameter and others that measure the temperature 
structure parameter. Some of the methods include analysis of 
star trails on photographic emulsions [Ref. 3], an 



2 



isoplanoxneter which measures the isoplanatic angle through 
stellar scintillation [Ref. 4], and a Modulation Transfer 
Function (MTF) device for determining the coherence length 
[Ref. N 5]. These devices measure important integrated 
parameters but they cannot measure the vertical profile of 
turbulence. An acoustic echosounder which is similar in design 
and construction to a SONAR [Ref. 6], measures a vertical 
profile but is usually limited to several hundred meters 
range. Greater vertical resolution occurs at the expense of 
decreased maximum range. 

In order to get a measurement of the vertical profile of 
turbulence to 20-50 km a device must be able to be launched 
on a sounding balloon or an aircraft. An example of this type 
of device is the thermosonde originally designed by GTE 
Sylvania and revised by the Air Force Weapons Laboratory and 
Tri-Con Associates Inc. and built and used by the Air Force 
Geophysics Laboratory [Ref. 7]. 

This thesis is an attempt to design and build a 
temperature sensing probe system to measure the vertical 
profile of the temperature structure parameter which is a 
measure of the vertical profile of turbulence. Although such 
systems exist, they are 1) expensive, greater than $2000 per 
launch and 2) require extensive calibration. The purpose of 
this thesis was to develop an inexpensive device (the system 
cost is approximately $150 per launch) which is simple to 
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operate (it is a self-calibrating device) . Additionally the 
effects of solar heating of the probes are investigated. 

The results of a comparison with a well developed 
acoustic echosounder indicates the system will be effective 
in measuring turbulence in the atmosphere. The studies of the 
solar heating effects indicate the only source of error due 
to solar heating would be from the sun/shade effect on the two 
probes and if that is corrected for, it will be accurate up 
to 30 km. 



4 



II 



BACKGROUND 



A. TURBULENCE PARAMETERS 

Small temperature variations carried by the turbulent 
velocity field in the atmosphere produce small phase 
perturbations in an optical plane wave propagating through it. 
These perturbations randomly distort and convolve the phase 
of a plane wave. There are three atmospheric parameters which 
must be determined prior to any attempt made at compensating 
for these atmospheric distortions. These parameters are the 
refractive turbulence structure parameter, C N 2 , the spatial 
coherence length of the atmosphere, r 0 , and the isoplanatic 
angle 0 O . 

The most important of these parameters is C N 2 . Tatarski 
[Ref. 3] states that one way to deal with a non-stationary 
problem, which includes all atmospheric parameters, is to 
define a function in terms of a difference, he than defines 
the structure function for index of refraction as, 

D n ( r ! , r 2) = < [ N( r 2 ) - N( r j ) ] 2 > / (1) 

where < > denotes an ensemble average. If we assume the 
atmosphere to be homogeneous and isotropic over small regions 
the structure function can be rewritten as. 
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( 2 ) 



D„(r) = <[N (r 2 ) -N (r^) ] 2 >, 

where r is r^-r z . By dimensional analysis Kolomogorov [Ref. 3] 
showed that the structure function has an r 2/3 dependency. 
Consequently D N is proportional to a constant C N 2 times r 2/3 , 

D n = C N 2 r 2/3 . (3) 

C N 2 is the refractive turbulence structure parameter, a mean- 
square statistical average of the difference in the index of 
refraction between two points separated by r 12 , 

C N 2 = <(N 2 -N 1 ) 2 >/r 12 2/3 . (4) 

The r 2/3 normalization extends from an inner scale 1 0 , on the 
order of millimeters, to an outer scale of L 0 , on the order of 
meters [Ref. 3]. These fluctuations in the index of refraction 
arise from variations in density caused by temperature 
fluctuations in the turbulent velocity field. 

These density variations in the atmosphere alter the 
phase of an optical beam being propagated through it. The 
Optical Transfer Function (OTF) characterizes the integrated 
phase perturbations of an optical beam. It is a measure of the 
correlation of the electric fields of the optical beam 
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perpendicular to the direction of propagation. Although the 
atmosphere is not homogeneous or isotropic, Tatarski [Ref. 8] 
postulates the idea of local homogeneity and isotropy, which 
states that over some region R, comparable to the outer scale 
length L 0 , we can assume the atmospheric random variables are 
homogeneous and isotropic. The modulus of the OTF is the 
atmospheric Modulation Transfer Function (MTF) . Fried [Ref. 
9] introduces the parameter r 0 to characterize the MTF. It is 
related to the refractive turbulence structure parameter by. 



r o = 21 



1.46k 



t C n 2 < z H 



-3/5 



(5) 



where r 0 is the spatial coherence length, k is the wave 
number (2*/A), and C N 2 is the refractive turbulence structure 
parameter along the optical path of length L [Ref. 5] . 

The other measure of spatial coherence in the atmosphere 
is the isoplanatic angle 0 O . It is similar to r 0 in that it 
is the dependence of the optical transfer function of a system 
for different angles to the source. The two parameters r 0 and 
6 0 are conjugate pairs. d 0 looking up is equivalent to r e 
divided by the path length looking down [Ref. 4] and vice 
versa. A more formal definition is that 8 0 is an angular 
measure of spatial coherence, it is the limiting angle for 
which an electromagnetic wave from a source will follow the 
same optical path length to a receiver. 
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Figure 1. Comparison of LASER Beam Propagating 

Through Vacuum and Atmosphere, showing r 0 
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Figure 2. Two Sources Showing the Effects of 
the Isoplanatic Angle. 
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If we consider two paths through the turbulence, the 
isoplanatic angle relates the mutual coherence e' 1 point of the 
E field between the two paths. Fried [Ref. 9] expresses the 
isoplanatic angle's dependence on C M 2 as, 



*o = 



2.905k" 



C n 2 (z)z*/3 dz 



-3/5 



( 6 ) 



where z is the altitude. Note the z 5/3 spherical weighting 
factor that emphasizes turbulence far from the optical system. 

The preceding paragraphs clearly show the importance of 
the C M 2 parameter, it not only defines the measure of 
turbulence, it also determines the spatial coherence r 0 and 
the isoplanatic angle 8 0 . The problem lies in that high 
resolution profiles of C N 2 are difficult to measure, it 
requires complex detectors and optical imaging systems. 
Instead we can define a temperature structure parameter C T 2 
similar to C N 2 where, 



C T 2 = <(T 2 -T 1 ) 2 >/r 12 2/3 , 



(7) 



which can be measured by several different methods. The 
fluctuations of the index of refraction are due to 
fluctuations in the density of the atmosphere and if we can 
assume that the density fluctuations are due solely to 
temperature fluctuations, then C T 2 is related to C H 2 by, 
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( 8 ) 



Gp 2 , 

where P is the atmospheric pressure in millibars and T is the 
atmospheric temperature in Kelvins [Ref. 10] . The assumption 
that index of refraction fluctuations are due only to 
temperature fluctuations and that humidity fluctuations are 
insignificant, is valid when the Bowen Ratio B (ratio of 
sensible heat flux to latent heat flux) is greater than 0.3 
[Refs. 11,12]. Below this value, humidity fluctuations are 
significant. 

B. SCALE LENGTHS 

In Kolomogorov 's definition of the structure function he 
assumed local homogeneity over a region bounded by the inner 
and outer scale lengths 1 0 and L 0 . These scale lengths vary 
from hundreds of meters at the outer scale lengths down to 
millimeters for the inner scale lengths [Ref. 13]. The outer 
scale length is the size of the turbulent fluctuations at the 
onset, while at the inner scale viscosity dissipates the 
energy of turbulence as heat. Kolomogorov called the region 
between 1 0 and L 0 the inertial subrange and as long as the 
distance r, in Equations (4) and (7), is within this inertial 
subrange these equations are valid. Kolomogorov expressed the 
power spectral density of the turbulence in this region by. 



c n 2 = 



79x10' 6 P 



T 2 
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( 9 ) 



<t>(K) = 0.033C 2 K' U /3 
n 



where 2*L 0 ' 1 << K « 2 jt1 0 ' 1 



Von Karman took this definition 
further to include the ranges for eddy sizes greater than L 0 , 



r( 1 1/6) x ' 9 / 2 

^ = "r( i/3) 8 



<<5 2 > L 3 
n o 



- 11/6 



1 + 



K 



V 



( 10 ) 



where K 0 =2*L 0 ' 1 and <5 n 2 > is the variance of the refractivity 
fluctuations and is related to C N 2 by. 



C H 2 - 1.9<6 n 2 >K 0 2/3 , 



(ID 



[Ref. 10]. Figure 3 shows the spectrum of turbulence for all 
scale lengths, the inertial subrange shows Kolomogorov' s 
linear description of turbulence while in the regions above 
2n 1 0 1 viscosity effects dominate and below 2*rL 0 ' 1 Von Karman's 
spectrum defines the turbulence [Ref. 14]. 

It is important to understand the effects of the 
turbulence spectrum in all three regions because of the errors 
introduced due to the lack of correlation in the regions above 
and below the inertial subrange. Care must be taken in 
choosing the correct value of r to cover the scale lengths 
inside the inertial subrange. Additionally the temporal 
frequency of the turbulent fluctuations is related to the wave 
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number K and the wind speed moving the turbulence past the 
probes. Therefore to measure the thin transition layers in the 
atmosphere accurately [Ref. 15] the errors due to various 
scale sizes must be determined. These calculations will be 
carried out in Chapter 4. 



SPECTRUM OF TURBULENCE 

FOR ALL SCALE LENGTHS 




Figure 3. Spectrum of Turbulence Over All Scale 

Lengths . 



C. SYSTEM REQUIREMENTS 

A fast response time temperature probe attached to a 
sounding balloon system can measure the vertical profile of 
the temperature structure parameter. The design of the 
temperature probe system provides a simple, low cost method 
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that is capable of resolving the thin stratified layers of a 
stable atmosphere and determining the temperature gradient, 
thickness, and turbulence of these layers. In designing the 
system the key points of consideration were 1) what type of 
probe geometry and thereby what type of circuitry and 2) what 
type of temperature sensing element to use. 

The probe system could use a differential measurement or 
a single point probe. The advantage of the single point probe 
is that it can measure C T 2 from either the variance of the 
data, knowing the balloon ascent rate, or by analyzing the 
power spectral density .However , this would require a high data 
rate (several samples per second) , which the radiosonde 
systems used do not have. Although a differential system would 
not have as high a vertical resolution it has the advantage 
of providing partially reduced data which relaxes the need for 
a high data rate. Therefore a differential system greatly 
simplifies the telemetry needed for the system at the expense 
of more complexity in the sensor itself. The resolution for the 
system would still be satisfactory, about 2 meters of vertical 
resolution for a balloon with an ascent rate of 2 meters per 
second . 

There are several different choices for the probes. They 
ican be made from resistance wires, thermistors, or 
thermocouples . The resistance wire is simply a fine wire, such 
as platinum or tungsten, with a known resistance as a function 
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of temperature. The problem with a resistance wire for this 
application is that it requires some means of self-calibration 
due to its dependence on Ohm's law and temperature. A 
thermistor is a small semiconductor device which changes its 
resistance as a function of temperature. It has a larger 
change in resistance vs. temperature then a metal wire, 
although it is non-linear. However, thermistors are large 
compared to a fine wire, having a larger thermal mass than a 
probe made from a fine wire, which increases the response time 
and the susceptibility to solar heating of the device. 
Additionally, both the thermistor and the resistance wire 
require a current source which not only increases the 
complexity of the circuit but also introduces a self heating 
factor. A thermocouple consists of two fine wires made of 
dissimilar metals welded together, which produces a voltage 
difference proportional to the temperature. For an in-depth 
discussion of thermocouples and the thermoelectric effect, see 
Refs. 16 and 17. Commercially available thermocouples can be 
made of very fine wire (down to 12.5 ^m) thereby reducing the 
thermal mass and producing a faster response time. A 
disadvantage of a thermocouple is that the response to 
temperature is small, typically 40 nV/°K. This places severe 
requirements on low noise signal processing. An operational 
requirement for this type of device would be the knowledge of 
the mean atmospheric temperature to calculate the Seebeck 
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Coefficient, which is the derivative of the thermal emf with 
respect to the temperature. Since the rawinsonde system used 
with the probe provides this data, it is easily accomplished. 
Small temperature changes of 1° C or less, are expected from 
the two probes in a differential system. This produces a 
negligible change in the Seebeck coefficient therefore no 
calibration of this type of device is needed, other than 
knowing the gains of the electronic amplifiers. 



\ 
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III. SYSTEM DESIGN AND DEVELOPMENT 



A. AMPLIFIER CIRCUIT 

The probe system consists of a pair of thermocouples 
connected in series and held rigidly 1 meter apart by an 
aluminum tube. An amplifier circuit requires high voltage gain 
and very low noise to be able to discriminate the signal from 
the background due to the low voltage produced by 
thermocouples (on the order of n volts per degree) . It also 
requires an analog root mean square device to calculate <(T 2 - 

Ti) 2 >. 

Since the purpose of the system is to measure temperature 
differences on the order of hundredths and thousandths of a 
degree and the thermocouple output voltage is about 40 
microvolts per degree C the circuit must have a very high gain 
with ultra low self noise. Therefore the circuit must be 
carefully designed and built to reduce any sources of noise 
wherever possible, by methods such as matching resistor values 
as exact as possible and placing them as close as possible to 
each other to reduce the thermal drift. Other examples include 
thermal insulation for the entire circuit and RF shielding for 
the circuit as well as the probe. The most critical component 
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of the circuit is the ultra low noise, precision, high speed 
op amps which have voltage noises less than those of 50 ohm 
resistors. 

The circuit (designed by Prof. Don Walters and fabricated 
by Dale Galarowicz) is a low noise, wide bandwidth 

Instrumentation amplifier. Three operational amplifiers 
produce a gain of 10,000 and a combination low and high pass 
filter with a gain of 5 producing a total voltage gain of 
50,000. The circuit contains two Linear Technologies LT1028 
ultra low noise precision high speed op amps in the 
instrumentation amplifier portion of the circuit. These op 
amps have a gain bandwidth product of 75 MHz and a self noise 
of 0.85 nV/Hz 1/2 at the frequencies desired. Figure 4 shows the 
noise and frequency characteristics of this op amp. The filter 
for the circuit uses an LT1057 op amp whose noise and 
frequency response characteristics are shown in Figure 5. 
Based on the frequency response curves in Figures 4 and 5 
[Ref. 18] and the high pass filter the circuit has a frequency 
response from 0.16 Hz to 200 Hz. The other major component of 
the circuit is an Analog Devices AD637 high precision, wide 
band RMS-DC converter. It has a bandwidth of 600KHz at lOOmV 
RMS and an averaging time constant of 25 msec/micro F. The 

entire circuit is powered by two 9 volt dry cells. Figure 6 

\ 

is a complete schematic diagram of the circuit. 
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Voltage Noise vs Frequency 
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Figure 4. Noise and Gain Characteristics of the 

LT1028 OP AMP. 
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Figure 5. Noise and Gain Characteristics of the 

LT1057 OP AMP. 
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Figure 6. 



Schematic Diagram of Differential 
Amplifier. 
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The circuit has to be well shielded against RF 
interference due to the high gain of the circuit and the fact 
it will be operated in the atmosphere where it is highly 
susceptible to all types of RF signals. Additionally the 
thermocouple wires must be stretched out over a distance of 
1 meter and will act as an antenna. The aluminum tube shields 
them from RFI except at the two endpoints, greatly reducing 
extraneous signals. The circuit is inclosed in a styrofoam 
casing, to reduce thermal gradients across the circuit, which 
is then covered in aluminum foil. The circuit itself has 
ferrite beads and an LC filter, at the inputs, to further 
shield the op amps from RFI and the entire circuit is built 
on a ground plane which has been grounded with the foil 
shield. 

The result of the aluminum foil and styrofoam shielding 
and use of ultra low noise op amps is an amplifier capable of 
measuring the extremely small voltages produced by the 
thermocouple probes. The circuit was tested in the anechoic 
chamber in the basement of Spanagel Hall to minimize any 
temperature fluctuations and then measurements of the circuits 
self noise were taken. Figure 7 shows the noise spectrum 
measured by a Hewlett Packard HP3561A Dynamic Signal Analyzer. 
The large noise spike below 1 Hz is due to 1/f or flicker 
noise which occurs in all amplifiers, due in large part to 
surface leakage of transistors [Ref. 19]. Although the noise 
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normally occurs at frequencies up to 100Hz it has been reduced 
by the use of a high pass filter. The noise spike at 60Hz is 
due to AC powered equipment operating in the anechoic chamber. 
This noise should vanish when the system is used in the field, 
since it is powered by DC batteries, as long as care is taken 
to insure it is not near a large AC power source. The plot 
clearly shows the self noise output of the circuit is well 
below 100 microVolts/Hz 1/2 in the frequencies of interest. 
Since the circuit gain is 50,000 and the Seebeck coefficient 
is on the order of 40 microvolts per degree this translates 
into a noise induced measurement of less than 0 . 00005°C/Hz 1/2 . 



RANGE: 1 dBV STATUS: PAUSED 

A: MATH SORT (MAG' 2 / BW) RMS: 100 




Figure 7. Noise Spectrum of Amplifier Circuit. 
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The entire system was set-up and run in the anechoic 
chamber with caps on the exposed thermocouples to further 
reduce the signal. After running for about 2 hours to let all 
air currents in the chamber settle, the system produced the 
results shown in Figure 8. This shows the noise produced by 
the circuit introduced an error corresponding to a C T 2 of less 
than 10' 6 , which is two orders of magnitude lower than the 
lowest C T 2 needed for a usable probe system. These results 
indicate the circuit self noise is well below that which would 
have a degrading effect on the results. 

The entire circuit package measures 3" X 3" X 3" and 
weighs less than six ounces. 

B . THERMOCOUPLE 

In 1821 Thomas Seebeck discovered that two wires of 
dissimilar metals joined at one end and heated produce a 
voltage difference across the open ends. This voltage is a 
function of the junction temperature and the composition of 
the two metals. Since that time, many different thermocouple 
types have been produced based on the combination of the two 
different metals used and having different thermoelectric and 
physical properties. The requirements for the system included 
use in the atmosphere from the surface to 20 km, therefore the 
temperatures vary from 30° C to -30° C with as large a Seebeck 
coefficient as possible. The Seebeck coefficient is the slope 
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of the voltage versus temperature curve at a given 
temperature. Based on these requirements the T type or Copper- 
Constantan thermocouple was selected. An alternative choice 
would be the E type or Chromel -Const antan thermocouple. 

The Copper-Constantan thermocouple is composed of a 
copper wire and a 55% copper 45% nickel wire. It has a 
temperature range of -200° C to 350° C and is suitable for 
applications where moisture is present. Table 1 shows the 
thermoelectric voltages referenced to 0° C for a Copper- 
Constantan thermocouple based on the National Bureau of 
Standards reference tables. 

9 

Table 1 shows that the Seebeck coefficient, which is the 
unit difference in voltage for each temperature change, is not 
linear over the entire temperature range. To determine the 
Seebeck coefficient, dV/dt, valid over the entire temperature 
range desired, the data over the temperature range desired 
from Table 1 was plotted and then a polynomial regression was 
performed to find the equation of the curve. The derivative 
of the curve was taken to find the Seebeck coefficient. A 
fifth order polynomial fits the data from -100° C to 30°C. The 
Seebeck Coefficient for a Copper-Constantan thermocouple is, 

dV/dt=3 . 8707x10 2 +8 . 5348xl0' 5 t-3 . 3135xl0‘ 7 t 2 

-2 . 77432x10 9 t 3 -l. 253xl0 11 t\ (12) 
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Table 1 



VOLTAGES FOR A TYPE T THERMOCOUPLE 
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where t is in °C and V is in millivolts. The temperature 
difference measured by the probe, is the output voltage 
divided by the circuit gain and divided by dV/dt determined 
at the specific temperature. For a 1° change in air 
temperature the average change in dV/dt is less than 0.2 
microvolts per degree C or less than 0.5% change in dV/dt. 



Cons tan tan 




Figure 10. Schematic of Thermocouple Probe. 

The probe consists of two 0.00254 cm diameter Copper- 
Constantan thermocouples held rigidly 1 meter apart by an 
aluminum tube. This wire has a resistance of 998.3 ohms per 
double meter and the circuit requires a low input impedance, 
therefore it is necessary to cut the fine thermocouple wire 
short and solder larger wires to them to run the distance 
between the thermocouples and the circuit input. Copper- 
Constantan wires 0.0254 cm in diameter were used which have 
a resistance of 9.983 ohms per double meter. It is essential 
to solder the copper to copper and constantan to constantan 
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to insure no other thermocouple junctions are formed which 
could cancel any voltage signal generated. The resulting probe 
yields an input resistance of 30 ohms. 



1/TRU VS SQRT of WIND SPEED 

FOR Cu-Cons t an t an THERMOCOUPLE .00254cm DIAMETER 




Figure 11. Response Time Versus Wind Speed of a 

.00254 cm Copper-Constantan Thermocouple . 

The reason for the fine wire was to reduce the thermal 
mass which decreases the response time of the thermocouple to 
temperature change. The response time is defined as the time 
required to reach 63.2% of an instantaneous temperature 
change. The OMEGA Handbook [Ref. 20] gives the response time 
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for a 0.00254 cm thermocouple as 0.05 sec in still air and 
0.004 sec in 18.3 m/sec air. The response time r can be 
related to the velocity by the equation, 

\ 

_1 = a + iy~y- , (i3, 

T 

where A and B are constants for the specific thermocouple 
material size and are determined by the response times given 
in Ref. 20 (for this case A=20 and B=53.783 sec' 1 ). This 

results in a response time verses velocity curve shown in 
Figure 11. From this curve the response time r § 5 m/sec is 
.0071 sec. The interaction of the frequency response and the 
inner scale size will be discussed in the Results section. 

C. COMPUTER AMD CODE 

The data acquisition and processing portion of the system 
consists of a Hewlett Packard model 217 computer with a 20 
megabyte hard drive and 2 megabytes of memory. It also 
contains an Infotek BC203 Basic compiler, an Infotek AD200 
analog to digital converter, an HP9133 floppy disk drive, a 
monitor and printer. The computer receives data from the 
circuit via the analog to digital converter. It reduces the 
data and displays a C T 2 verses time (altitude) profile. 

The analog to digital converter is the primary component 
for the data acquisition, however it introduces errors. The 
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A to D converter has a small DC voltage offset. This offset 
is small so it does not significantly affect the results until 
C T 2 reaches values on the order of 10' 3 , as is discussed in the 
results section. Each A to D converter has a different offset 
so each time a different computer is used the offset must be 
measured. In addition the noise of the electronic circuit 
produces a DC offset. To find the total DC offset, the entire 
system was set up and run in the anechoic chamber to insure 
no external signals are introduced. The software has a feature 
for inputting the offset and another to display the voltages 
directly from the probe. To find the offset, initially set the 
offset to zero and then after the system has run for about one 
hour, to settle all the air currents generated by the set-up, 
read the voltages by pressing the PRINT_RAW key, this displays 
the actual voltage offset then simply reboot the system and 
input the offset read from the raw voltages. The DC offset for 
the system used for the experiment was -5 millivolts. 

Appendix A summarizes the features of the software and 
lists the code for controlling and producing output from the 
system. It is almost fully automatic. Once it has started 
running, the only corrections to be made are to update the 
ambient air temperature any time there is a change of two or 
three degrees, so the Seebeck coefficient can be recalculated. 
There is a hard function key for updating the air temperature. 
A further refinement to the system would be to have the 
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rawinsonde temperature measurement device automatically update 
the program with the new air temperature at a given interval 
to make the system fully automatic. 



D. SOLAR HEATING 

It is important to understand the effects of thermal 
radiation on the probes since they will be ascending through 
the atmosphere both during the day and at night. During the 
day radiative heating from direct and reflected solar 
radiation will heat the probes above ambient air temperature 
while at night the probes will be cooled due to Planck 
radiation from the probes to space. There are several ways in 
which the heating or cooling of the probes can introduce 
errors into a differential measurement device. A hot wire 
anemometer effect can introduce errors. Another is the 
difference caused by one probe being in the sun and the other 
being in the shade or at night where one probe has a direct 
line of sight to the earth and the other is blocked as by a 
cloud. In a hot wire anemometer effect the probes are heated 
above the ambient temperature and then velocity fluctuations 
across the probes vary the heat transfer rates away from the 
probes, creating a false temperature difference [Ref. 21]. In 
the sun-shade effect, as the probe ascends through the 
atmosphere one end may come into shade either from a cloud or 
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the balloon shadow. The probe may rotate as it ascends giving 
one thermocouple a different aspect to the sun than the other, 
creating the same sun/shade effect to a lesser extent. In 
either case one thermocouple will receive less direct solar 
radiation than the other thereby changing the heat flow on one 
thermocouple but not the other, again introducing a false 
difference. The same applies at night if one thermocouple is 

shaded from the earth by a cloud it will not radiate thermal 

energy at the same rate as the other will. 

The net heat flow to or from a body in the atmosphere is 
described by the heat transfer equation, 

q ne t=q s + <3a + <J t + 3c + - <3r * ( 14 ) 

where q s = portion of direct solar radiation absorbed 

q a = portion of the atmospheric radiation absorbed 
q t = portion of terrestrial radiation absorbed 
q r = thermal radiation emitted by the wire 
q k = net conduction to the wire from the atmosphere 

q c = net convection to the wire from the atmosphere 

[Ref. 2.2]. The temperature difference can be determined by 
first setting q net equal to 0. The equation can then be reduced 
into the heating and cooling portions by, 

^SD + ^SR + ^LE = ^lu + E c , (15) 
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where E^, = heating due to direct solar radiation 

E sr = heating due to solar radiative reflection of 
v the atmosphere 

E le = heating due to long wave radiation from earth 
E Lg = cooling due to long wave radiation from wire 
E c = convective cooling 

[Ref. 23]. If we further assume the thermocouple to be a 
horizontally oriented, infinite cylinder with the top half 
radiating to the sky and the bottom half radiating to the 
earth's surface, the temperature difference between the 
ambient air and the thermocouple is given by, 



AT = 



j a 



1 



R s + Tf L 



R + 

a g ^4 

o 



(16) 



where h = average convective conductance 

e s = short wave emissivity of the thermocouple 
£ L = long wave emissivity of the thermocouple 
a = albedo 

R $ = short wave incoming radiation 

R fl = long wave radiation from the earth's surface 

R a = long wave atmospheric radiation 

o = Stephan-Boltzmann constant 

T = air temperature 

where h is defined by Kreith [Ref. 22] as. 
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h = Kk 

D - 

where D = wire diameter 
v V = wind speed 

k = heat conductivity of the air 
v = kinematic viscosity of the air 
K&n = empirically determined dimensionless constants 
based on the Reynolds number 

[Ref. 24]. 

The Air Force Geophysics Laboratory [Ref. 23] defines the 
calculation for delta T in a similar manner based on the same 
assumptions. However, when comparing the two forms using the 
same parameters (see Figure 12) there is a considerable 
difference. Review of both treatments shows the principle 
difference lies in each definition of the convective 
conductance h. Campbell [Ref. 24] uses Krieth's form of h (EQN 
17) while Brown [Ref. 23] assumes an average value. Another 
method of determining h is with the Nusselt number, which is 
a dimensionless number used in describing heat transfer and 
fluid flows. Kramers [Ref. 25] performed extensive 
measurements of heat transfer on spheres and cylinders, from 
this he determined the Nusselt number to be a function of the 
Reynolds number and the Prandtl number, another dimensionless 
number where, 



VD 

v 



( 17 ) 
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Re = VD p / p. , 



( 18 ) 



and Pr = c^/k , 

where D = diameter of the cylinder 
V = velocity of the fluid 
p = density of the fluid 
p = dynamic viscosity of the fluid 
k = thermal conductivity of the fluid 
c p = specific heat of the fluid. 

From all the available data he showed that the Nusselt number 
for a cylinder could be represented by, 

Nu = 0 . 91 (Pr) 0,31 (Re) 0-385 , 0.1<Re<50 , 

and (19) 

Nu = 0.60(Pr)°‘ 31 (Re) 0 ' 50 , 50<Re<10 , 000 . 

Based on the Nusselt number the convective conductance is, 

h = Nu*k/D. (20) 

Campbell included some experimental results in his paper. 
When the experimental results are compared with Kreith's and 
Kramers' treatment of h we can see that Kramers' expression 
exactly models the actual data (Figure 13) . Using Equation 20 
for h in Equation 16 corrects the differences seen in Figure 
12, therefore the lower curve in Figure 12 is the correct 



34 



model for solar radiative heating in the atmosphere, based on 
the results obtained using the experimental data of Reference 
24 . 

SOLAR HEATING OF THERMOCOUPLE 



o 



_J 

u 




Figure 12. Comparison of AFGL Model and Campbell's 
Model of Solar Heating of Thermocouple 
Wires . 



A discrepancy noted in Campbell's calculations was the 

values used for short (visible) and long(IR) wave emissivities . 

\ 

Table 2 gives the emissivities used and the actual 
emissivities from a 1986 edition of the CRC handbook. The 
corrected values were used to recalculate the solar heating 
and the updated results are shown in Figure 14 . 
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Figure 13. Comparison of Convective Conductance. 



Table 2 

EMISSIVITIES OF CU-CONSTANTAN 
AND TUNGSTEN 





Visible 


IR 


Campbell's 

Cu-Constantan 


.25 


.5 


CRC Handbook 
Cu-Constantan 


.2 


. 03 


CRC Handbook 
Tungsten 


.5 


. 03 
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SOLRR HERTING OF THERMOCOUPLE 
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RLTITUDE (Km) 

Figure 14. Corrected Comparison of solar heating. 



Based on Figure 14 it is now possible to make a 
determination of the errors introduced by solar heating. The 
structure function for velocity fluctuations over small scale 
lengths is, 



N D (r) = 3 . 83 (« r) 2/3 , (21) 

where e is the dissipation rate (Ref. 26]. Actual data shown 
in Reference 26 from areas of highly turbulent velocity fields 
indicates the average dissipation rate is on the order of 
3x10 5 m 2 sec' 3 . For the probe system with r equal to 1 meter 
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this yields velocity fluctuations on the order of 0.06 m/sec. 
Figure 15 shows the solar heating errors due to the hot wire 
anemometer effect with velocity fluctuations of this 
magnitude. Figure 16 shows the differences between the curves 
of Figure 15. It indicates the temperature differences are 
negligible (on the order of .001 degrees or less). 

The shading effect can be determined by eliminating the 
direct solar radiation component in the equation. This will 
show the maximum error, if it is a matter of the probe 
changing aspects to the sun the errors will be proportionally 
less. If direct solar radiation is completely removed there 
will be very slight heating of the probe, due to incoming 
terrestrial radiation. The temperature difference between the 
two probes will be approximately equal to the amount of solar 
heating on one probe as seen in Figure 14. This indicates a 
major source of error since the temperature difference at 
higher altitudes is on the order of 0.2 °C and the balloon can 
rotate as it ascends. 

Based on these calculations if some method is devised to 
limit the rotation of the probe to eliminate the sun/shade 
effect, the probe system can measure C T 2 values accurately up 
to 30 km altitude without significant errors. 

Reference 7 described a C T 2 thermosonde used by the Air 
Force Geophysics Laboratory. Data measured by this system 
indicates an order of magnitude jump in the values of C N 2 
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DELTR T -(C) • DELTA 



SOLAR HEATING OF THERMOCOUPLE 




ALTITUDE (Km) 

Figure 15. Hot Wire Anemometer Effect on Copper 
Constantan Thermocouples. 



DIFFERENCE BETWEEN THE CURVES 




ALTITUDE (km) 

Figure 16. Differences Between the Curves of Figure 

15. 



39 



just after sunrise [Ref. 27]. This increase appears to be an 
artifact of the instrument rather than actual turbulent 
processes due to the fact it occurs so rapidly. Since it 
occurs at sunrise a logical assumption is that it is due to 
solar heating, therefore a great deal of time has been spent 
in determining these effects. Figure 17 shows the hot wire 
anemometer effect on a 4 micron tungsten resistance wire. The 
values for the emissivity of Tungsten are taken from Table 2 . 
This indicates solar heating does not effect the measurements 
of Cj 2 consequently the rise in the value of C n 2 at sunrise may 
be actual. 



SOLAR HEATING OF 4 MICRON TUNGSTEN WIRE 




Figure 17. Hot Wire Anemometer Effect on 4 micron 

Tungsten Wire. 
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IV. RESULTS 



A. EXPERIMENTAL PROCEDURE 

The experimental measurements served two purposes, first 
they were carried out to validate the probe measurement system 
and second, they were used to validate the C T 2 measurement 
capabilities of the acoustic echosounder [Refs. 6,28]. The 
acoustic echosounder calculates a 15 minute time averaged 
value for the temperature structure parameter as a function 
of altitude, however uncertainties in the antenna beam shape, 
side lobes, transducer efficiencies and atmospheric 
attenuation produce uncertainties in the absolute value of C T 2 
calculated by the echosounder [Ref. 28]. Therefore independent 
verification of the C T 2 values must be made to validate the 
acoustic echosounder. If both values agree this is a positive 
indication that both systems are measuring accurately. 

The acoustic echosounder is a high frequency device which 
analyzes the atmospheric density fluctuations within the first 
200 meters of the atmosphere. The echosounder consists of a 
Hewlett Packard HP 217 computer to control the system and 
acquire and reduce the data, an HP 3314A function generator 
which produces the pulsed signals, an amplifier and the array 
of speakers which acts as a transmitter/receiver. The system 
operates at 5KHz and produces a 100 cycle sinusoidal burst of 
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18 acoustic watts. The antenna array consists of 19 
piezoelectric speakers in a close-packed hexagonal array 
enclosed in a 55 gallon plastic trash container lined with 
lead cind foam to reduce side lobes as much as possible. The 
minimum range of the device is approximately six meters based 
on the recovery or "ring" time of the speakers and the maximum 
range is about 200 meters based on the frequency used. Figure 
18 is a diagram of the acoustic echosounder layout. 

A comparison test was run with the probe system and the 
echosounder between 1300 and 2030 hours local time on 3 
September 1988 on the roof of Spanagel Hall on the grounds of 
the Naval Postgraduate School at Monterey CA. The acoustic 
echosounder was placed on the sixth level, northwest corner 
of the roof while the probe was attached to a rigid pole and 
extended approximately 1.5 meters off the seventh level of the 
northwest corner of the roof. In this position the probe was 
approximately 9.2 meters above the echosounder array, thereby 
being outside of the echosounder blind zone. Due to the 
building itself and heating exhaust vents on the eastern side 
of the building it was necessary to monitor the wind direction 
to insure the prevailing winds were not passing over the 
building and picking up heat from the exhausts, which would 
have greatly affected the data. Therefore wind speed and 
direction as well as temperature and humidity measurements 
were taken every 15 minutes to update temperature and humidity 
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FOR FILTER, USE 
INPUT GAIN I ANO 
OUTPUT GAIN 10. 



Figure 18. Layout of Acoustic Echosounder 
Device. [Ref . 28] 
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information for the two systems as well as determining if the 
prevailing winds were flowing across the building before 
passing over the instruments, corrupting the data. During the 
entire experiment the wind shifted several times but it was 
always from Northwest to Southwest. At all times the air flow 
passed over the instruments before passing over the building. 
All the data was valid. 

Figures 19 through 21 represent a portion of the data 
collected during this experiment. In each of the figures the 
upper graph is the echosounder profile of the atmosphere, the 
central graph is a 15 minute time averaged C T 2 profile of the 
atmosphere based on the data collected by the acoustic 
echosounder and the bottom graph is the C T 2 versus time plot 
measured by the probe system at an altitude of 9.2 meters 
above the echosounder. The dark lines below about 6 meters in 
the upper plot and the discontinuities below 6 meters in the 
center plot are due to the blind zone of the echosounder. To 
compare the plots, the average of the bottom plot was compared 
with the value at an altitude of 9.2 meters in the center 
plot. Figure 19 was taken early in the afternoon and shows 
strong convective pluming, which causes a higher temperature 
structure parameter. Figure 20 which was taken closer to the 
neutral event shows a marked decrease in the turbulence and 
a corresponding decrease in the values for C T 2 . Figure 21 
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Figure 19. 



Echosounder Trace and C T 2 Measurement and 
Probe Measurement During Strong Turbulence. 
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Figure 20. Echosounder Trace and C T 2 Measurement and 
Probe Measurement During Light Turbulence. 
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Figure 21. 
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TIME AVERAGED CT'2 





Echosounder Trace and C T 2 Measurement and 
Probe Measurement During No Turbulence. 
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taken during the neutral event, shows virtually no turbulence 
and a much lower value for C T 2 . Appendix D contains additional 
measurements taken during this experiment to show the 
corresponding increases and decreases in C T 2 for both the 
echosounder and the probe system. The purpose of these 
measurements was not to actually sample the atmospheric 
processes at this location but to make a quantitative 
comparison test between the acoustic echosounder and the 
probe system. For a complete description of the atmospheric 
turbulence measurements and processes for this location see 
Weingartner [Ref. 6]. 

B. SCALE SIZE ERRORS 

The temperature variations in a turbulent atmosphere range 
in size from millimeters to hundreds of meters. Optical 
aberrations are primarily caused by variations the size of a 
Fresnel zone (AD) 1/2 therefore with laser frequencies and path 
lengths of several kilometers the important scale sizes are 
on the order of several centimeters . [Ref . 13] With a frequency 
response of 150 Hz and an average wind speed of 2-5 m/sec, the 
system is limited to scale sizes greater than 3 cm, which will 
introduce a small amount of "inner scale" error, since minimum 
scale sizes are on the order of millimeters, but if used for 
measurements in conjunction with laser propagation through the 
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atmosphere the error will be negligible. Additionally, the 
acoustic echosounder utilizes the smaller scale sizes of 
approximately 3 cm, thus there will be negligible error 
introduced by this in a comparison test. 

To find the outer scale length errors we can express the 
structure function of the probe system with limiting scale 
lengths by, 



D(a,b) = 4 



(wrj) 2 s i n 2 ( k r / 2 ) k“ 5 / 3 



dk 



a l + ( wr j ) 2 l + ( ut 2 ) 2 



( 22 ) 



where a is the limiting lower frequency (outer scale) 

b is the limiting upper frequency (inner scale) 

k = 2II/A where A is the actual scale length 

w = kV where V is the wind velocity 

t , = RC time constant of the high pass filter 
(the upper frequency cutoff) 

r 2 = frequency response of the probes 
(the lower frequency cutoff) 

r = probe separation distance 

and' comparing it with the structure function over all scale 
lengths, 



D(0,oc) = 4 



XX) 




. 2 


' k r ' 


sin 




- 0 


~2 





-5/3 



dk , 



(23) 



the outer scale length limiting error can be determined. 
Figure 22 is a graphical representation of this comparison, 
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showing the error over a range of limiting outer scale 
lengths. The 8% limiting error on the low end of Figure 22 is 
due to the high frequency cutoff of the circuit and as the 
scale size decreases the larger low frequency errors of the 
circuit begin to dominate, increasing the error. The limiting 
scale length for the experiment can be determined as the 
height above ground, which was approximately 30 meters. From 
Figure 22 it is clear the error introduced due to finite outer 
scale lengths is approximately 12%. The design of the acoustic 
echosounder is resistent to outer scale errors therefore the 
finite inner and outer scale error introduced would cause the 
probe system to record measurements approximately 12% lower 
than the echosounder. 
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The acoustic echosounder is susceptible to inner scale 
errors which will cause it to read higher than the probe 
system. The inner scale length is inversely proportional to 
the wind speed and can be expressed as, 

1 0 = 7.4(„ 3 /0 V4 . (24) 

Ochs and Hill [Ref. 29) made extensive measurements of the 
inner scale length, based on their results and the mean wind 
speed of 6 m/sec during the measurements, the approximate 
inner scale length was 3 mm. At the edge of the inner scale 
of turbulence there is a bump in the temperature spectrum due 
to diffusion as it enters the viscous-convective range. Figure 
23 illustrates this bump showing the spatial power spectrum 
4 > t of temperature fluctuations versus the scaled wave number 
k f? , which is the wave number normalized by the inner scale 
length. [Ref . 30] Here k is equal to 2w/Scale Length and r? is 
equal to 1^7.14 (for air). The limiting inner scale size of 
the echosounder is 3.4 cm (A/2 where A = (340m/sec)-i- (5kHz) ) 
therefore with an inner scale size of 3mm the scaled wave 
number is approximately 9xl0’ 2 . Figure 23 shows the acoustic 
echosounder will read approximately 5% higher than the 
Kolomogorov spectrum and therefore 5% higher than the probe 
system. 
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K’l 



Figure 23. Spatial power spectrum $ T of temperature 
fluctuations versus scaled wave number ktj . Solid curve 
is actual model; the dashed curve is Tatarski's model. 
[Ref. 30] 



C. ANALYSIS OF DATA 

Figure 24 is a comparison of 15 minute time averaged data 
collected from the echosounder and the probe system. This data 
was taken before the noise measurements and discovery of the 
5 millivolt DC offset error in the A to D converter RMS module 
combination. Figure 25 shows the corrected data comparing the 
two systems. The data clearly shows the correlation of the two 
systems even with the volatile trends of the turbulent 
fluctuations. It also shows a decrease in the temperature 
structure parameter leading up to and during the neutral 
event, which corresponds with the actual physical processes 
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Figure 24. Comparison of Data From Acoustic 

Echosounder and Temperature Probe Before 
Correction for DC Offset. 
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Figure 25. Comparison of Data From Acoustic 

Echosounder and Temperature Probe After 
Correction for DC Offset 
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TABLE 3 



C T 2 MEASUREMENTS 
(CORRECTED FOR DC OFFSET) 



TIME (LOCAL) 


ECHOSOUNDER 


PROBE 


% DIFFERENCE 


1400 


. 0354184 


. 0368390 


3 


1415 


. 0316086 


.0391012 


23 


1430 


.0358514 


. 0335893 


6 


1445 


. 0500970 


.0492531 


2 


1500 


.0322137 


. 0321643 


0 


1515 


.0490100 


. 0432243 


12 


1530 


.0474621 


. 0384319 


20 


1645 


.0311984 


.0271304 


13 


1700 


.0263566 


.0350252 


32 


1715 


.0361568 


. 0292298 


20 


1730 


.0260007 


. 0232678 


10 


1745 


. 0217524 


. 0274032 


25 


1800 


. 0233976 


.0200511 


14 


1815 


. 0141386 


. 0165680 


17 


1830 


. 0138539 


.0165680 


18 


1845 


. 0132466 


. 0184468 


39 


1900 


.0129296 


.0145711 


15 


2000 


. 0054777 


. 0077345 


40 


2015 


. 0058262 


.0071837 


23 




AVERAGE % DIFFERENCE 


17 



going on at the time. With the approach of sunset, at 1933 
local time, the sun heated the earth's surface to a lesser 
degree thereby reducing the temperature difference between the 
earth's surface and the air, which in turn reduced the 
temperature structure parameter and the turbulence. 
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Table 3 contains the values of C T 2 measured by each of the 
devices and corrected for the offset, it indicates an average 
difference of 17%, with the probe system reading lower. The 
scale length errors indicate the probe should read 
approximately 12% lower due to outer scale length errors and 
5% lower due to inner scale bump errors. Therefore there is 
no significant difference between the probe system readings 
and the echosounder calculations. 



56 



V. CONCLUSIONS 



Independent verification of C T 2 values measured by the 
acoustic echosounder is important [Ref. 28] and the 
differential temperature structure parameter probe has 
provided a valuable comparison indicating the absolute C T 2 
values of both the echosounder and the probe are valid. Taking 
into account all known errors there is no significant 
difference between readings of the echosounder and the probe 
system, which is an extremely good indication that both 
systems are providing valid measurements. Additionally this 
thesis demonstrated that solar heating of the probes in the 
atmosphere does not appear to play as significant a role as 
first thought. The only major effect solar heating has on the 
differential system is when one probe is directly illuminated 
by solar radiation while the other is shaded. 

This probe system has many applications including being 
placed on towers to calibrate other turbulence measuring 
devices as well as being attached to a rawinsonde system and 
launched to measure the vertical profile of C T 2 . When used in 
this mode it can measure values of C T 2 up to 30 km altitude 
accurately, however if it rotates as it ascends through the 
atmosphere the sun/shade effect of solar heating will 
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adversely affect the system. Further developments to the 
package, such as addition of wind vanes on the probe assembly 
which will not affect the turbulent flow but will dampen the 
rotation, or a small motor with a flywheel to act as a 
gyroscopic stabilizer to prevent the probe from rotating, will 
eliminate the errors induced by this effect. 

Other improvements to the system include methods to 
automatically update the temperature into the program from 
the balloon systems onboard temperature sensor. Another 
improvement would be to increase the data transmission rate 
of the rawinsonde telemetry system to get a higher resolution 
profile of the thin stratified layers of the turbulent 
atmosphere or even possibly having the system transmit an AC 
signal from which a great deal more information can be 
extracted such as the power spectral density of the 
turbulence. 
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APPENDIX A 



PROBE SYSTEM SOFTWARE 

The program that runs the system is called "CTSQR". It 
controls the probe system, collects the data, reduces it, and 
then displays and stores it for further analysis. It is based 
on the same program that controls the acoustic echosounder. 
It can be broken down into several sections. The first section 
sets up the system, initializes all arrays, creates a data 
file which can store up to eight hours of data and sets up the 
function keys which are used to update the temperature used 
in calculating the Seebeck Coefficient, prints out raw voltage 
data or ends the program storing what has been collected. The 
next section initializes the Infotek AD200 analog-to-digital 
converter, which collects the data. Now that the system is 
ready to collect data it calculates the Seebeck Coefficient 
based on the information input at start-up or updated through 
the function key. Next it collects data every second and 
averages it over ten seconds, reduces it to C T 2 and plots it 
every ten seconds. Every 15 minutes it prints out the plot and 
then resumes the data collection. The program "CTREADER" can 
take the data file generated by "CTSQR" and read it and 
calculate 15 minute time averaged values of C T 2 . 
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IF Npo int <Npo 1 nt_ol d THEN Store_dete«l 




930 


IF Npoi nt_ol d>0 AND Store_dete-0 THEN 




940 


1 Plot the dete 




950 


Ctsqr«( Vol t ege/t Ge 1 n+Seebeck «R_one_ 


.third) )*2’ 


960 


Lgct sqr-L6T { Ct sqr ) 




970 


* MOVE Npolnt_old,Volts_old 




900 


■ DRAW Npolnt , Vol tage I Plot 


Vo 1 1 eoe 


990 


MOVE Npolnt_old,Lgct eqr_ol d 




100e 


DRAW Npolnt ,Lgc t eqr 




1010 


END IF 




1020 


Vcl t s_ol d-Vol t ege 




1030 


Lgct sqr_ol d-L get sqr 




1040 


Npoint_ol d-Npolnt 




1050 


IF Store_dete-l THEN Npolnt_old-0 




1060 


NEXT 1 




1070 


Volt s_avg-SQR( Vol tsq/10 ) 




10e0 


Ct-Vol ts_evg/( 6ein«Seebeck*R_cne_thlrd ) 




1090 


01SP ‘ ‘ 




1 100 


T1-T1ME0ATE 




1110 


Deyf-OATEt ( T 1 ) 




1 120 


Time$-T1ME»(T1 ) 




1130 


i 




1140 


i ••• CT SQUARED DATA REDUCTION SECTION 




1 150 


1 




1160 


Yrf-DayJie, 11 1 




1 170 


Dey-( (PATE<Dey*)-OATE< "1 JAN *&Y r* ) ) 01V 06400)+! 


1190 


1 




1190 


Nrec-Nrec+ 1 




1200 


« CALCULATE DECIMAL HOURS 




1210 


T?*T inaf 




1220 


Hcur--VALt Tft 1 , 2 )) + < VAL( T$t 4 , 5 1 )+VAL( Ttl 7 , 8 1 )/60 )/60 


1220 


Hr* 1 000*Hcurs * Note that the HP rounds 




1240 


PRINT 'Record t’lNrec;" Collected *;Deyt; 


’ ’ ; T 1 ne$ 


1250 


PRINT ‘ ‘ 




1260 


ALPHA OFF 




1270 


• SET UP OUTPUT ARRAY 




12e0 


Nl-Nrec+1 IFOR OPTION BASE 1 




1290 


02< 1 , 1 )*Nrec 




1300 


02<N1 , 1 )-Oay 




1310 


02 ( fil .2 >-Hr 




1320 


02( N1 , 3 )*Vo 1 1 e_avg* See 1 e 




1330 


0?( N 1 ,4 )-C t • See 1 e 1 Intensity 




1240 


PRINT N 1 ; 02 ( N 1 , 1 ) ; 02 ( N 1 , 2 ) ; 02 ( M 1 , 3 ) ; 02 ( N 1 


,4 ) 


1350 


Store_riata'. < Write output every NPLDT seconds 


1260 


IF Stcre_dete-1 THEN 




1370 


S t cre_det e-0 




13P0 


PRINT 




1390 


01SP 'WRITING P.FOUCE 0 OUTPUT’ 




1400 


OUTPUT PFi lei ;D2< « ) 




1410 


ASSIGN PFilel TO Filelt 




1420 


P 1 O t nun*P 1 o t nun + 1 




1 4’0 


IF P 1 ct nun MOO 2-1 THEN 




1440 


PRINTER IS 701 




1450 


PRINT ’ 
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1460 

1470 

nep 
1< 90 
15P0 
1 51 P 
152P 
1530 
1540 

isse 

15E0 
1570 
1590 
1590 
1G00 
1610 
1G20 
1G30 
1B4P 
IE50 
1660 
1670 
16 eP 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
I7E0 
1770 
1790 
1 79P 
iep0 
1910 
ie2P 
ie:p 
1P40 
1950 
ie60 
1 97P 
1990 
ie90 
1900 
1910 
1920 
1930 
1940 
1950 
I960 
197P 
1900 
1990 
2PPP 
2P10 
2P2P 
2P3P 
2P4P 
2P5P 
2P60 
207P 
2PPP 
2P90 
2 1P0 
2110 
2120 
2130 
2140 
2 15P 
2160 
2170 
2190 
2190 
2200 
22 IP 
222P 



PRINTER IS 1 
ENO IF 

DUMP 6RAPHICS *701 1 Dump screen to printer 

Npc i nt -0 
DISP ' ' 

CALL Plot setup( Nplol , Si ta* ) 

END IF 

End_whi I a END WHILE 

I 

L*pdat e_t enp ■ < Updates Seeback Coefficient With Maw Air Temp 

INPUT 'ENTER NEW AIR TEMPERATURE( DEGREES C)“,T 
See“3 . 9707E-2+9 . 534BE-5*T-3 . 31 35E-7»T*2 
Beck— 2.77432E-9*T'3-1 .253E-1MTH 
Seabeck-( Sea + 9eck ) • 1 . E — 3 
GOTO Star t_lo 

Prlnt.raw. I TD6GLE THE PRINT FLA6 
Pri nt_key— Pr i nt_key 
GOTO Start_lo 

Quit- * Write reduced data output flla 
FOR 1-1 TO Nl 

PRINT 1:02(1,1 ) ;02< 1 , 2 ) ; D2< I,3);02< 1 ,4 );D2M ,5);02< 1 ,6 ) 
NEXT I 

OUTPUT PFllel ;D2 < • ) 

PRINT 'DATA FILE HAS PEEN STORED UNOER NAME " , F 1 1 a 1 * 

BEEP 

BEEP 

ASSIGN OF 1 1 a I TO * 

STOP 

END 

SUB PI ct se tup( Np lot , Si t eS , Ymln ) 

Ymax** I 
GINIT 

6RAPHICS ON 

LINE TYPE 1 

VIEWPORT 15, 120, IS. e0 

WINDOW 0 , Np 1 ot , Ym l n , Ymax 

AYES 60, .5,0, Ymln, 5, 2 

CLIP OFF 

CSIZE 4, .6 

LOPG 6 

{ Draw Log Y A:< 1 s 
1 FOR Oecade-Ymln TO Ymax 
i FOR L‘nlts-1 TO MB* < Oacade<Ymax > 

' Y-Decade+LGT (Units ) 

i MOVE 0 , Y 

« DRAW Nplot.Y 

t NEXT Units 

« NEXT Decade 
1 Label horizontal axis 

TMTIMEDATE MOD 96400 
Hrs-Tl DIV 3600 
T2-T1 MOO 3600 
Ml n“T2 DIV 60 
Qtrhr»Mln DIV 15 
FOR M*0 TO Np 1 ot STEP 300 
MOVE M, Ymln-. 06 
Qtrmln-Qtrhr*15+( M/300) *5 
IF Otrmln-60 THEN 
Ot rml n-0 
Hr s-Hr s4 1 
ENO IF 

LABEL USING “00 , A , Z 2 * ; Hr s : “ • * :Qtrnln 
NEXT M 

MOVE Np I ot /2 , Yml n- . 0 
LAeEL 'Tima (Local )' 

1 Label Ordinate 
LORG 9 

FDR M-Ymln TO Ymax 
LORG 8 
CSIZE 4 

MOVE -Nplot/23 . 3 , M 
LABEL USIN6 *t,K";"10“ 

CSIZE 2 
LORG 1 

MOVE -Mplct/2e,M 
LABEL USING ‘*,K*;M 
NEXT M 
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rt'bt? 

2000 
2100 
2110 
2120 
2130 
2140 
2150 
21G0 
2170 
21B0 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2250 
2270 
22P0 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2350 
2370 
23P0 
2390 
2400 
2410 
2420 
2430 
24 40 
2450 
2450 
2470 
2400 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2550 
2570 
25P0 
2590 
2500 
2510 
2520 
2530 
2540 
2550 
2550 
2570 

2Ge0 

2590 

2700 

2710 

2720 



fit X I fl 

MOVE NpIot/2,Ynln-.D 
LAPEL ‘Tina (Local >" 

1 Label Ordlnata 
L0R6 B 

FOR M-Ynln TO Ynax 
LORG B 
CS1ZE 4 

MOVE -Npl ot /23 . 3 , M 
LABEL USIN6 , » 1 K*: , 10‘ 

CS1ZE 2 
L0R6 1 

MOVE -MpJct/2B.M 
LAPEL USIN6 *»,K“;M 
NEXT M 
LOIR Pl/2 
L0R6 5 
CS1ZE 4 

MOVE -Nplot /7, < Ynln+Ywas )/2 
LABEL "L06 OF CT*2" 

(Title the p 1 ot 
LOIR 0 
LORG 4 

MOVE Hplct/2,Ynax+l 

LABEL OATES ( T 1ME0ATE >; “ ";Sltet 

CLIP ON 

SUBENO 

SUB Adi n( Vcl t age , Zero ) 

( 25 APR 1 9B5 OLU 

t INFOTEK A-0 Input rout lna set up for lntarnal trigger 
• and average 40 points ovar three 50 Hz cycles 
! 

INTEGER 1 , Npoint a, Ad_dat a( 1 -.40) 

P1M S e 1 ec t S[ 20 I 
Ad_;el_ccde-17 
Int ensl ty*0 
Npolnts-40 

Count f-VALf < Npoint s ) 

Scale-5. 0/<Npolnts«2047. ) 

9 t de v-0 

De 1 t a_t lnef*" 1 250000" t AD Interval between aanples In naac 

Se 1 net * •" ae I act Isl end" 

G05VE Read_ad 

FOR 1-1 TO Npolnts 

Volt age-Vol t age *Ad_dat a( 1 ) 

NEXT 1 

Vo 1 1 age-( Vol t age- Zero )*Scal a 
GOTO Subend 

Read_ad ( Read the Infotak A-0 

• Initialize the A-0 

OUTPUT Ad_sal_coda; "RESET" , "internal " , "count ", Counts 
OUTPUT Ad_ael_code : " 1 1 ne " ,0a 1 1 a_t Inel, ‘dal ay on" , Sal acts 
OUTPUT Ad_seI_code: "STATUS" 

ENTER Ad_sel_code : Resp* 

IF RespS-" " THEN 

ENTER Ad_sel _code USIN6 " t,U" ; Ad_data< • ) 

OUTPUT Ad_aal_code; ‘STATUS- 
ENTER Ad_aaI_code : Resp* 

IF ResptO " " THEN 

PRINT "ERROR 0UR1NG SAMPLIN6 - ";Resp* 

ENO IF 
ELSE 

PRINT ‘ERROR 0URIN6 A-0 INITIALIZATION - *:Ra»pl 
ENO IF 
RETURN 
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2730 Subend- 1 

2740 SUBENO 

2750 SUB Init_ad200 

2780 • INITIALIZE AD_200 

2770 Code-17 

2780 Ourmy-RE AOIO< Code , 3 ) 

2790 WRITEIO Code, 0.0 

2800 CONTROL Code,0i I 

2810 SU8EN0 
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APPENDIX B 



SOLAR HEATING PROGRAM 



13 t 

23 

23 

43 

53 

S3 

73 

93 

93 

133 

113 

123 

123 

143 

150 

1S3 

173 

123 

123 

230 

213 

223 

223 

243 

2sa* 



RE-STORE 'ALTPLOT!.- ,730,1,3" 

• 12 OCT 1929 MRC 

* THIS PROGRAM PLOTS THE CHANGE IN THE DELTA T ' TEMPERATURE CIFFEPENCE 

* BETWEEN THE THERMOCOUPLE ANC THE AIR TEMP CUE TO SOLAR HEATING) 

• WITH INCREASE IN ALTITUDE 2ASE0 ON SCCC’S CALCULATIONS IN AF SL PL 9 
! CA.TEO 22 FES 1294 AND CAMPBELL'S WORK FROM CCT 12S9 



i THE FOLLOWING IS THE LIST OF VARIABLES 
! 

J-4.18 ’ 

S -.14 1 

Epot-.25 
Spaa-. 25 
Fd-PI 

Fr-2. ! 

Epaiw-.5 ' 

Epol-.S » 

Ra-. 322 ' 

Ra-.33S \ 

Rg-.315 1 

Albedo-. 25 * 

S i gma~5 . S7E- 12 » 

Sigcai-1 . 2SS39E- 
Dn*. 3333254 ■ 

Den-. 30254 i 



nechanicai equivalent : f heat (W Cal“-1 See) 

Solar Constant (W cn*-2 ) 

Viaible Wavelength ABSORPTIVITY of the THERMOCOUPLE 

Short Wave Emaaivity fer CAMPBELL'S Calculat icna 

Fern Factor for Cirect Radiation 

Fern Factor for Reflected Radiation 

Lsng Wave EMISSIVITY ef tha THESflCOCUPLE 

Long Wave Eni aa: v : iy for CAMPBELL’S Calculations 

Short Wave Inoonirg Radiation (Cal cn*-2 aec*- 1 > 

Long Wave Atnoapheric Radiation (Cal cn A -2 aeo*-i > 

Long Wave Radiation from Ground (Cal cn*-2 aec*-l) 

Reflection ef the Garth 

St efan-Bo 1 1 :wrn Constant (W cn*-2 K*-4> 

12 J Stefan-8ol t rnann Ccnatant (Cal cn A -2 K*-4) 

Probe Dianeter in netera 
Probe Dianeter in centinetera 
)E X 2 1 Denver a 1 or Factor for Thernal Conductivity 
( Conatant for Deternining Mu (Kq*- 1 n*-l K * — 1 / 2 > 

1 Sutherland'a Conatant for Mu (K> 

( Olnenaionleaa Prandtl Nunber for Air 
Cinenaioni eaa Reynolda Nunber for air 
Dinenai on leaa Nuaaelt Nunber for A:r 
Air Temperature (Kelvtna) 

Air Denaity aa Function of Altitude (kg/n*2> 

Dynamic Viacoaity of Air (N-aec/n‘2) 

Kinenatic Viacoaity of Air (,V2/aec> 

Thernal Conductivity (Cal aec*-l cm*-2 (C/cn>*-l> 
Convective Heat Tranafer Coefficient (f(ALT)} 

(Cal cn*-2 aec"-l K‘-l > 

Percentage of Solar Radiation Reaching Given Altitude 
Heat Flux fron Earth at Given Altitude (W cm*-2> 

The 6iven Altitude (Kn> 

Wind Speed (n/aec> 

Enpirical conatant baaed on Reynolda Nunber fron Krieth 
Empirical conatant baaed on Reynolda Nunber fron Krieth 
Tenoerature Difference between TC and Air for GOOD 
Temperature Difference between TC and Air for CAMPBELL 
Both in Degree a C 

K> SPEED CD YOU WANT - ,V 



N PLOTTER IS 707, # HPSL" 

512 \ BE SURE PLOTTER DIP SWITCHES ARE PROPERLY SET ie. SWITCH 1,2,2 IN PCSIT 1 
514 IF 0-2 THEN PLOTTER IS CRT , * INTERNAL * 

523 > SINIT 

523 » GRAPHICS ON 

S40 LINE TYPE 1 

553 VIEWPORT 15,120,15,90 



2S3 


Ccnver 


273 


Beta-1 


283 


Suth-1 


223 


Pr-. 71 


233 


i Re - 


213 


'Nud - 


223 


)T 


323 


^Rcw - 


343 


\ Mu - 


253 


»Nu - 


263 


* K 


373 


»H 


393 




223 


»F 


439 


'Se - 


419 


lAlt - 


420 


1 M m 


423 


1 Kr ~ 


449 


IN 


459 


! Del t - 


469 


! Del 1 1- 


470 


> 


480 




490 


INPUT 


- S00 


l 


513 


t PLOT 


SI1 


INPUT 


512 


IF 0-1 
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552 'J INCOW 0,32,0,1 

570 AXES 5,. 05,0 0,5,2 

520 CLIP OFF 

580 C5I2E 4 .5 

£00 LCRS S 

SI0 » LAS EL HORIZONTAL AXIS 

520 FOR M-0 TO 30 STEP S 

520 MOVE M, - . 04 

£40 LA2EL M 

550 NEXT M 

£50 MOVE 15,-. i 

570 L.ASEL "ALTITUDE (Km>" 

5 80 i L.ASEL VERTICAL AXIS 

590 LCRS 9 

702 FOR M-0 TO 1 STEP . 1 

710 MOVE - . 5 , M 

722 LA2EL M 

720 NEXT M 

740 LOIR PI/2 

752 LORE S 

750 MOVE -4.5 .5 

770 LABEL "DELTA T CC>" 

780 * TITLE PLOT 

790 LOIR 0 r 

800 LORS 4 

919 MOVE 15,1.1 

929 LASEL "SOLAR HEATINS OF THERMOCOUPLE" 

939 MOVE 15,1 

940 C5IZE 2 

959 LASEL "FOR A U I NO SPEED OF " ; V ; " m/ aee " 

950 CSIZE 4 

972 !• 

999 I- Calouiat :cna 

999 ! 

909 Aiteid-0 

912 Celtcid-0 

929 Oeltleld-0 

929 FOR Alt-0 TO 29 

940 \ I atm aaauning the rel.it tcnahip for F end Se ere linear urt Altitude 

950 f The Values were taken from eROWN and 6CCD 

950 F-. 5+1 A 1 1 * .912903 > 

970 Se-. 056-1 Alt* .00141935 > 

999 \ The relat icnabipa fcr T ,Rcw, Mu.Nu , and K are taken from the 

999 * HAND9CCK of GEOPHYSICS and the SPACE ENVIRONMENT chap 14 

1000 !■ and baaed on a U.5. Standard Atncaphere 

1010 * Thia ia taken from the Standard Atmosphere Temperature Profile 

1022 IF Alt <-10 . THEN 

1939 T-289.15-?.015*Alt 

1040 ENO IF 

1050 IF AltM0. AND Ait<29 THEN 

1059 T-21S . 

1970 END IF 

1099 IF Alt >-20 . THEN 

1090 T-219.+1.2MAit-20. > 

1120 END IF 

1119 Rcw-1 .2252-. 0399*A1 1 ! Denaity Change ia approximately Linear wrt 

1129 * Altitude up to about 100 Km 

1130 Mu-(8eta*T A 1.5 ).MT+Suthl 

1140 Nti-f Mu/Rcu >*. 1224255 1 Convereicn faotor to get correct viacoaity 

1150 Power— 12. /T 

1150 K1-2.£52195-2*T“1.5 

1170 K2-T+245 . 4 *19 "Power 

1190 K-<K«/K2 1 /Convert 

1192 f The determination of H waa baaed on KRAMERS (194S) which waa ahown 

1209 \ tc mere clcaeiy approximate experimental reauita by the plot HPL0T2 

1212 Re~Rcw*V*On/Mu 



66 



1223 IF Ro-<50 THEM 

1220 Nud m . 21 ^Pr* .31 •R^'* .225 

1240 EL5E 

1250 Nud~.S*Pr A .2l*Ro“.5 

1250 ENO IF 
1272 H-Nud^K/Cc* 

1223 i Colculatisns basod :n SRDWN and 5O0D ( AF5L Pub 1224? 

1252 1 

1200 ' PRINT "ALT " ; A1 i , “ROW* ; Rc^ f "MU* i flu t "NU" ; Nu 
1210 ' PRINT " K":K t "R&":Rs* t "H"!H 1 ,, T*iT 

*223 i PRINT “KRIETH’S H-“ ; Cl , “FIRSTPART-" ;A1 , “ 2ND-" ; ei 
1320 A-l./iJ*HJ 

1240 9-F*S*Ep:>tM ( 1 ./Fd >+( Albsds/Fr ) } 

1250 C-EpsIwM (Se/Fri-lSisn^T**) > 

12E0 Dol t~A*< 2+C 1 r 

1270 LINE TYPE 9 

1220 MOVE A 1 1 o 1 d , Del t s I d 

1220 CRAW A lt,Dsli 

1402 ! Calcuiaiisno based on CAMPBELL U9S9) 

1410 l 

1420 IF Re<4 THEN 

1420 Kr-.SSl 

1440 N-.220 

1450 END IF 

14S0 IF Re<40 AND Re>'4 THEN 

1470 Kr*^. 822 

1490 N m . 2 2 5 

1420 END IF 

1503 IF Re>-40 THEN 

1510 Kr-.SIS 

1520 N-.46S 

1520 _END IF 

1542 7tl-£pss*Rs*< l + < ( P I *Albeds 1/2 >) 

1550 21-PI*EpalM URa+Ra>/2 >-<Sia::ai*T A 4 ) > 

1553 Cl-<Ds:n/{Kr*K > >MNu/<y*Drr>> A N * Krieth’s empirical for* of 1/K 

1570 Deltl-<A1+51 )+C 1 

1520 LIME TYPE 4 

1520 MOVE AI is id.Del tlcl d 

1500 DRAW Alt .Celt 1 

1510 Altold-Alt 

1520 Del t lold-Oelil 

1520 DeltcId~DeIt 

1540 PRINTER 15 701 

1550 PR INT "ALT'“ : Alt f “CELT ( AF5L ; Co-1 i , “DELT( CAMPBELL s Del i 1 

1550 NEXT Alt 

1570 D£ I ZE 2 

1520 LINE TYPE 9 

1522 MD'JE 1 .24 

1720 DRAW 2 , . 94 

1710 DRAW 2, .94 

1723 DRAW 4,. 24 

1720 MOVE 7 .94 

1740 LINE TYPE 1 

1750 LA2EL “AFGL FORM" 

1750 LINE TYPE 4 
1770 MOVE 12 .94 
172e DRAW 21 , .94 
1720 LINE TYPE 1 
1200 MOVE 25 ,.24 
1210 LAS EL “CAMPBELL FCRM“ 

1220 DL’MP GRAPHICS 4701 
1223 END 
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APPENDIX C 



2 » 

a r* 

50 

c r \ 

u U 

70 

20 

OA 



1 ^ n* 

i. *- u 

130 

140 

i cr* 

X ~JXJ 

1 cr* 
i <j v 

i 72 

1 20 

1 O'* 

A 

200 
n i a 

211 
-» i 

n i 4 
lie • 

220 ■ 
TTO» 1 

4. _/ XJ 

2 40 
250 
2 S 0 

‘-»*7£ 

290 

220 

300 

T i m 



340 

TCJ 

360 

370 

320 

3 90 

A (%f* 
-txjxj 

A 1 r* 

T A XJ 

420 

nnt 
-r v> 

440 

450 

4C« 
•» w w< 

470 

4 20 
4 90 
500 



CONDUCTANCE PROGRAM 

oc-cmqc -uoi ny-j. -mm t r% •* 

• ' k_ w « w n h_ n« w i - * ' wvi i i u 

« THIS program compares the results of the computation of the convective 

' CONDUCTANCE BY G. CAMPBELL’S METHOD AND BY THE AFGL PAPER’S METHOD 
! AND PLOTS THEM AGAINST EXPERIMENTAL VALUES FROM CAMPBELL 

I ii nry i ooo 

i VALUES ARE FOR A STD ATMOSPHERE la. 1 ATM AT 20 DEG C 

i 

3c;«j -1 . 25 1 k 5 / n " 3 AIR DENSITY 

-1 . 71 E-5 ! < N-sac >/:n"2 DYNAMIC VISCOSITY OF AIR 

n*2/5 KINEMATIC VISCOSITY OF AIR 
THERMAL CONDUCTIVITY OF AIR 
n DIAMETER OF THE PROBE 

DIMENSIONLESS PRANDTL NUMBER OF AIR AT 20 DEG C 
KK ANO N ARE EMPIRACAL NUMBERS BASED ON THE REYNOLDS NUMBER 
AS DESCRIBED IN KRIETH 19S5 pg 411 
V IS THE WIND SPEED IN r,/a 
Ro- IS THE DIMENSIONLESS REYNOLDS NUMBER 

Nud IS THE DIMENSIONLESS NUSSELT NUMBER AS DESCRIBEO IN KRAMERS 194S 

» CORRECTION FOR THE EMPIRICAL CALC BETWEEN KRAMERS 



Nu -1 . 722E-S 
K-S. 175E-5 
0--2.54E-2 
Pr 714 
! 



:rr3Ct". ,'bJ 



*7 C 1411P11C 



I 



• PLOT SETUP r 

INPUT “INPUT 1 FOR PLOTTER OR 2 FOR CRT* , Q 
IF Q-l THEN PLOTTER IS 707, "HP GL" 

IF Q -2 THEN PLOTTER IS CRT INTERNAL" 

BE SURE PLOTTER DIP SWITCHES ARE PROPERLY SET :a. SWITCH 1,2,3 IN POSIT 1 
GINIT 

GRAPHICS ON 
LINE TYPE 1 

iiTCi.ipnoT tc i *■>'* ic o a 

V i W \ x I A W y 1 C U f 4 W y W O 

WINDOW 0, 10,0, .5 

AVCC 1 1 r* f* C 

i t . * , t; , t; , j 

CLIP OFF 
OS I Z E 4 . S 
LORG S 

1 LABEL HORIZONTAL AXIS 

cno M-n> to i r* 

I Wit I < "<U 1 V-/ l XJ 

MOVE M ~ 02 
LABEL M 
NEXT M 
MOVE 5, -.07 
LABEL “VELOCITY Csn/sJ" 

i label vertical axis 

LORG 9 

COD M-m TO C CTCD I 

■ un i • XJ i \j . w i i— i . * 

MOI (C - A M 
i . w v t , i , 

LABEL M 
NEXT M 
LOIR PI/ 2 
LORG c 

MOVE - 1 . 5 , .25 

LABEL "h Teal cr."-2 see''-! C A -13* 

1 T T T L c PLOT 

i n t d rh 

UW 1 M XJ 

i noc a 
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510 


MHMC C _ C 


C~>n 


LA9EL "COMPARISION OF 


520 


MOVE 5 . 4 


540 


CSIZE 2 


CCn i 

*->«-> t-j 


LA8EL "WITH CORRECTION 


550 


PCT 7 C ^ 


570 


1 


con» 

wl W XJ 


He 1 d -0 


590 s 


1 i„ i rj-n 

V W A \J ~ XJ 


500 


MOVE 0 0 


c i ™ 

u 1 V 


LINE TYPE 5 


c * 

tl -T V 


FOR V“1 TO 10 


£50 


\ / m . . 

IvC'llUVD 1 V 1 1 U( A XJXJ • l t 1 •»-* 


CCA 

uuu 


TC TLlCM 

A I nGv^WWW 1 1 lull 


C rro 

«-l 1 XJ 


K “ 0229 


C Of* 

uuU 


m — one 

1 i • WUJ 


c 90 


GOTO Calc 


700 


END IF 


710 


re Ra> 4000 THEN* 


720 


K k - . 174 


720 


N - . 5 1 9 


n a n 

1 T U 


GOTO Calc 


750 


ENO T F 


7 C 0 


IF P. a N 4 0 THEM 


770 


Kk £15 


790 


. 4 £ £ 


-7Q n 
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